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Abstract. We study the effect of the core radius of an isothermal ellipsoidal gravitational lens on the caustics 
and the critical curves. We derive an analytic expression of the caustics for an isothermal ellipsoidal gravitational 
lens via a sixth-order algebraic equation. Since the expression is too long, by using another method we obtain a 
parametric representation of the critical curves in order to show analytically that there exist three cases: There 
are two curves for a small core radius, one for a quite large one, and no curves appear for an extremely large one, 
though the latter two cases are not realistic. The caustics are represented also by the same parameter. 
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1. Introduction 

Gravitational lensing due to a galaxy is indispensable to 
probe its mass distribution and determine cosmological 
parameters (for reviews, Schneider et al. 1992). Cored 
isothermal ellipsoids are often used for modeling a galac- 
tic lens. Although such a model is quite simple, it enables 
us to understand many physical properties of the galac- 
tic lens. In addition, it fits with mass profiles implied by 
observations (e.g., Binney & Tremaine 1987). 

The caustics are curves on the source plane perpendic- 
ular to the line of sight. They are important to understand 
the gravitational lensing, particularly when estimating its 
probability. Furthermore, the number of images changes 
if a source crosses the caustics. Similarly, critical curves 
located on the lens plane play a crucial role; giant arcs are 
observed near these critical curves. 

Analytic expressions for the caustics and critical curves 
were found for a singular isothermal ellipsoidal lens 
(Asada et al. 2003) during an algebraic study of the grav- 
itational lens (Asada 2002, Asada et al. 2002). However, a 
cored isothermal ellipsoidal lens, which is more important 
and realistic than a singular one, seemed beyond an ana- 
lytic approach, because the equation for a cored isother- 
mal ellipsoidal lens is of much higher order and thus more 
complicated from an algebraic point of view (Asada et al. 
2003). The existence of a core makes a crucial difference 
even if the core is small. For instance, the number of im- 
ages for a singular isothermal ellipsoid is two or four, while 



it is one, three or five for a cored isothermal one (Asada 
et al. 2003). The main purpose of the present paper is to 
study the effect of the core radius on the caustics and the 
critical curves. 

This paper is organized as follows. In section 2, the lens 
equation for a cored isothermal ellipsoid is given in the 
form of a single sixth-order algebraic equation. The dis- 
criminant for this algebraic equation is also discussed. A 
parametric representation of the critical curves and caus- 
tics are obtained in section 3. Section 4 provides our con- 
clusion. 

2. Lens equation for an isothermal ellipsoidal lens 

We consider an isothermal ellipsoidal lens with ellipticity 
< e < 1/5 and the angular core radius c. The surface 
mass density projected onto the lens plane must be non- 
negative everywhere. This puts a constraint as e < 1. A 
tighter constraint e < 1/5 comes from the requirement 
that the density contour must be convex, which is reason- 
able for an isolated relaxed system. The lens equation is 
expressed as 
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^ii-e)ei + ii + e)ei + c^ 

where (3 — (/3i,/32) and 6 — {61,62) denote the positions 
of the source and images, respectively. 

For simplicity, we introduce variables as x = \/l — e6i , 
y = vT+T6'2, a = VI ~ e/^i and b = vT+T/32, so that 
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the lens equation can be rewritten as 
a — X \ 1 



and 



6 = 2/1- 



1 + e 



(3) 



(4) 
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First, we consider Eq. I^J), which is rewritten as 



iR + c^f + c\l-e') 



For off-axis sources, this lens equation is reduced to a 
single sixth-order algebraic one as (Asada et al. 2003) 
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The other component is uniquely determined by 

(1 - <i)hx 



(1 
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Since the right-hand side of this inequality is positive for 
e < 1, this is recast in a polynomial form so that we can 
find the allowed region for R in an algebraic manner. This 
polynomial expression is 

F{R) - ' '•2^2 ^ 
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> 0. 
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It should be noted that the positions of the images can be The discriminant D4 of the fourth-order polynomial F{R) 
found analytically for a source on a symmetry axis (Asada is obtained (e.g. van der Waerden 1966) as 



et al. 2003). However, it is sufficient to consider a case of 
off-axis sources to obtain expressions of the caustics and 
critical curves because these expressions are continuous 
functions of the source position. 

By the discriminant for polynomial equations, we can This is non-positive, so that F{R) = has two real roots 



= -27c8(l + e)4 
X f (1 - 3e -I- 3e2 



(13) 



determine a location where the number of real solutions 
changes (e.g., van der Waerden 1966), which corresponds 



as i? = (1 - ef - and c^/^{l + efl^ - c^. We must find 
the smaller root satisfying a non-negativity of R. In short. 



to a change in the number of images in the context of we compare between 0, (1-6)^-0^ and c'*/^(l-|-e)2/^-c2. 
the lensing. The discriminant for the sixth-order equa- 
tion is computed for instance by using the software 
Mathematica (Wolfram 2000). The result is a lengthy 



There exist four cases for Eq. I|12(l . depending on c as 

(1 - e)3/2 



for c < 



polynomial of 656 terms, which is not tractable. Hence 
we employ another method to obtain a parametric rep- 
resentation of the caustics and critical curves in a much 
shorter and thus more practical form. 

3. Critical curves and caustics 

The critical curves are expressed as the vanishing of the 
Jacobian of the lens mapping, namely 9(a, h)/ d{x, y) = 0. 
It becomes 



(7) 



{x^ +y^ + cy - [x'' + y' + 2C-' + e{x^ - y^) 



for 



for 



for 



(l-he)i/2' 
0<i?<c*/'(l + e)2/3 
{1-ef 
(1 - e)3/2 



< R, 



< c < 1 
2 



x(a;^ + y2 + c^)V^ + c^(l_e^) = 0. 

This is rewritten as 

{R + - (i? + 2c^){R + c^l^ + c2(l - e^) 
= ei?(i?-hc2)i/2 cos 20, 



(l + e)i/2 
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< R, 



(14) 
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where we introduced the polar coordinates {x,y) = 
(r cos (j), r sin 0) and defined R = . 

We will find allowed regions among r, c and e, which 
is due to -1 < cos20 < 1 in Eq. ©. In short, Eq. © 
implies 



Next, we consider Eq. (|10|l . which implies 
G{R) = ((i? + c2)2 + c2(l-e2))' 

-(i?(l + e) + 2c2)V + c^) 



< 0, 



(20) 



(i? + -{R + 2c2)(i? + c2)i/2 + c2(l - e2) 
< eRiR + c^y/'^, 



whose discriminant D'^ is obtained (e.g. van der Waerden 
1966) as 

D'^ = -27c*^(l-e)4 

2 



(9) 
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(21) 



Asada: Isothermal Ellipsoidal Gravitational Lens 



3 



This is non-positive, so that G{R) = has two real roots 
as i? = (1 + e)^ — and c^^^\l — ep/"^ — c^. We compare 
between 0, (1 + ef - and c^/^jl - ep/^ _ c^. We find 
three cases for Eq. (I20II as 



for c < 1 — e, 

for 1 — e < c < 1 + e, 

0<i?< {l + ef -c^, 

for 1 + e < c, 

no allowed region for i?. 



< i? < (1 + e)2 - c^ (22) 

(23) 
(24) 

Bringing together Eqs. (jTU-lUni) and Eqs. 
we find four cases for Eqs. © and 1)10(1 : 
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(1 - e)3/2 

7^ < c < 1 - e, 

(l + e)i/2 

c4/3|l-e|2/3 

c4/3(l + e)2/3 
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<i?< (l + e)2-c2 
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These cases show that there are two critical curves for 
a small core radius c < 1 — e, one for the quite large case 
1 — e<c<l + e, and they disappear for an extremely 
large one, 1 + e < c. This anomalous behavior where the 
critical curves disappear is physically understood as fol- 
lows. As the core radius becomes much larger, the lens 
system approaches an approximately constant mass sheet. 
Around c ~ 1 + e, the size of the core is comparable to 
the Einstein ring radius. As a result, such a lens system 
cannot produce multiple images, which means that there 
are no critical curves. Furthermore, the latter two cases 
of large cores are not realistic. Indeed, c ~ 1 implies that 
the core radius is comparable to the Einstein ring radius, 
which becomes about lOkpc for a galactic lens with 10^^ 
solar mass at a high redshift. The core radius of a galaxy 
is usually much smaller than lOkpc. 

It is convenient to consider Eq. ^ as an r-parameter 
representation of the critical curves, because we must solve 
Eq. (jHl with respect to r to obtain a 0-parameter repre- 
sentation. This r-representation becomes 



{x,y) = (±ry ,±ry — ^), 

where we defined 

h^[eR)-^ ( {R + c^f'^-{R + 2c^) 

+c\l-e'){R + ^)-^'^). 



(31) 
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Fig. 1. Caustics for a cored isothermal ellipsoidal lens 
for e = 0.1 and c = 0.2. Sources are located at (0,0.1), 
(0,0.5) and (0,0.7), denoted by the triangle, filled disk 
and square, respectively. 



The allowed regions for r = ^/R are given by Eqs. ((25|l - 
(I30|l depending on c and e. Substituting this representa- 
tion into Eqs. Q and Q, we obtain a representation of 
the caustics as 



air) = ±r^ 



1 



b{r) = ±r 1 - 



1 - e 

1 + e 
Vr2 + c2 



1 + h 
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(33) 
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As an illustration, the caustics and critical curves are 
shown for e = 0.1 and c = 0.2 in Figs. 1 and 2, respectively. 
We can see two curves in both figures as discussed above. 
The inner caustics in Fig. 1 look like an asteroid. However, 
it is not exact but approximate, since Eqs. H33(l and H34|l 
are not the equation for asteroids (Asada et al. 2003). 

4. Conclusion 

We have studied the effect of the core radius of an isother- 
mal ellipsoid on the caustics and critical curves. For this 
purpose, we have obtained a simple parametric represen- 
tation of the critical curves. The caustics have been ex- 
pressed also with the same parameter. It turns out that a 
nonvanishing core radius changes the shape and number 
of the caustics and critical curves. These curves are two 
for a small core radius c < 1 — e, one for a quite large case 
1 — e<c< 1 + e, and disappear for an extremely large 
one 1 + e < c, where c is the angular core radius and e 
is the ellipticity of the isothermal ellipsoid. However, the 
latter two cases are not realistic in astronomical situations 
as discussed above. 

We must perform numerical computations for various 
values of the parameters in order to obtain the caustics 
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Fig. 2. Critical curves for a cored isothermal ellipsoidal 
lens for e = 0.1 and c = 0.2. The images correspond to 
the sources in Fig. 1. 



and critical curves. The present expression is analytic. 

Hence, it enables us to save time and reach more accurate 
results, for instance in rapid and accurate estimations of 
a lensing event rate, fittings to observational data, etc. 
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